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Abstract: We investigate the associated strangeness production reactions 
j~^ ' pp -^ pAK~^ and pp — *■ pH'^K'^ within an effective Lagrangian model. The 

04 , initial interaction between the two nucleons is modeled by the exchange of 

TT, p, CO, and a mesons and the strangeness production proceeds via exci- 
tations of 7V*(1650), 7V*(1710), and 7V*(1720) baryonic resonance states. 
The parameters of the model at the nucleon-nucleon-meson vertices are 
Q"^ , determined by fitting the elastic nucleon-nucleon scattering with an effec- 

CN| ' tive interaction based on the exchange of these four mesons, while those 

>P , at the resonance vertices are calculated from the known decay widths of 

the resonances and from the vector meson dominance model. Experimen- 
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^ij , tal data taken recently by the COSY-11 collaboration are described well 

by this approach. The one-pion-exchange diagram dominates the produc- 
tion process at both higher and lower beam energies. The excitation of 
the A''* (1650) resonance contributes predominantly to both the production 
Oh, channels at near threshold energies. Our model with final state interaction 

effects among the outgoing particles included within the Watson-Migdal 
approximation, is able to explain the observed beam energy dependence of 
the ratio of the total cross sections of these two reactions. 
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1 Introduction 

In recent years, there has been a considerable amount of interest in the 
study of the associated strangeness production reactions in proton-proton 
(pp) collisions. This is expected to provide information on the manifestation 
of quantum chromodynamics (QCD) in the non-perturbative regime of en- 
ergies larger than those of the low energy pion physics where the low energy 
theorem and partial conservation of axial current (PCAC) constraints pro- 
vide a useful insight into the relevant physics ^. The strangeness quantum 

^Talk presented in the second symposium on threshold meson production in pp and pd 
interactions, Jagellonian University, Cracow, Poland, May 31- June 3, 2004. 



number introduced by this reaction leads to new degrees of freedom into this 
domain which are expected to probe the admixture of ss quark pairs in the 
nucleon wave function [^ and also the hyperon-nucleon and hyperon-strange 
meson interactions 15 • 

The elementary nucleon-nucleon-strange meson production cross sections 
are the most important ingredients in the transport model studies of the 
i^"'"-meson production in the nucleus-nucleus collisions, which provide in- 
formation on not only the initial collision dynamics but also the nuclear 
equation of state at high density El 13 El El HIB El 021 • Furthermore, the 
enhancement in the strangeness production has been proposed as a signature 
for the formation of the quark-gluon plasma in high energy nucleus-nucleus 
collisions ;13. ■14; . 

The measurements performed in late 1960s and 1970s provided the data 
on the total cross sections for the associated hyperon (y)-kaon production 
at beam momenta larger than 2.80 GeV/c (these cross sections are listed in 
Ref. il5|). With the advent of the high-duty proton-synchrotron (COSY) at 
the Forschungszentrum, Jiilich, it has become possible to perform systematic 
studies of the associated strangeness production at beam momenta very close 
to the reaction threshold (see, e.g., Ref. J16j for a comprehensive review). At 
the near threshold beam energies, the final state interaction (FSI) effects 
among the outgoing particles are significant. Therefore, the new set of data 
are expected to probe also the hyperon-nucleon and hyperon-strange meson 
interactions. 

A very interesting result of the studies performed by the COSY-11 col- 
laboration is that the ratio (R) of the total cross sections for the pp -^ pAK^ 
and pp — > pTi^K~^ reactions (to be referred as AK^ and J^^K^ reactions, 
respectively) at the same excess energy (defined as e = y/s — mp — mY — rriK, 
with nip, niY, and mx being the masses of proton, hyperon, and kaon re- 
spectively and s the invariant mass of the collision), is about 28_g for e < 
13 MeV [17^. This result is very intriguing because at higher beam energies 
(e ~ 1000 MeV) this ratio is only around 2.5. 

Several calculations have been reported ^1^1120] to explain this result. 
Assuming that the tt- and K- exchange processes are the only mechanism 
leading to the strangeness production, the authors of Ref. ^H] show within a 
(non-relativistic) distorted wave Born approximation (DWBA) model that 
while the AK^ reaction is dominated by the ivT-exchange only, both K- 
and TT-exchange processes play an important role in the case of T,^K^ re- 
action. Therefore, if the amplitudes corresponding to the two exchanges in 
the latter case interfere destructively, the production of S° is suppressed as 
compared to that of A. It should however, be noted that in Ref. |19j . K- 



and IT- exchange amplitudes are reported to be of similar magnitudes for 
both AK^ and Ti^K'^ reactions. 

Nevertheless, a conclusive evidence in support of the relative signs of 
IT- and K- exchange amplitudes being opposite to each other is still lack- 
ing. Furthermore, other mechanisms like excitation, propagation, and decay 
of intermediate baryonic resonances which play (see, e.g., |21l 1221 \2'6\ ) an 
important role in the strangeness production, have not been considered by 
these authors. In the calculations reported in Ref. j2Uj also the relative sign 
of K— and it— exchange terms is chosen solely by the criteria of repro- 
ducing the experimental data, although in this work the theory has been 
applied to describe a wider range of data (which includes the polarization 
transfer results of the DISTO experiment [2H and the missing mass distri- 
bution obtained in the inclusive K~^ production measurements performed at 
SATURNE |2H] apart from the ratio R). 

We have investigated the AK~^ and T,^K~^ reactions at near threshold 
as well as higher beam energies in the framework of an effective Lagrangian 
approach (EL A) j21l fl2. 28, 29] . In this theory, the initial interaction be- 
tween two incoming nucleons is modeled by an effective Lagrangian which 
is based on the exchange of the vr-, p-, uj-, and a- mesons. The coupling 
constants at the nucleon-nucleon-meson vertices are determined by directly 
fitting the T-matrices of the nucleon-nucleon {NN) scattering in the rel- 
evant energy region. The ELA uses the pseudovector (PV) coupling for 
the nucleon-nucleon-pion vertex which is consistent with the chiral sym- 
metry requirement of the quantum chromo dynamics PT]. In contrast to 
some earlier calculations jSEl; both {AK^ and T,^K~^) reactions proceed 
via excitation of the A^*(1650), A^*(1710), and iV*(1720) intermediate bary- 
onic resonance states. The interference terms between the amplitudes of 
various resonances are retained. To describe the near threshold data, the 
FSI effects in the final channel are included within the framework of the 
Watson-Migdal theory j27l I29j . ELA has been used to describe rather suc- 
cessfully the pp — > ppn^, pp -^ pmT~^ [221 EHj, PP -^ pK~^Y [22 [22| as well 
as pp -^ ppe~^e~ [HO] reactions. 



2 Description of the Model 

We consider the tree- level structure (Fig. 1) of the amplitudes for the as- 
sociated K^Y production in proton-proton collisions, which proceeds via 
the excitation of the A^*(1650)(i"), iV*(1710)(i"^), and iV*(1720)(|"^) in- 
termediate resonances. To evaluate these amplitudes within the effective 



Lagrangian approach, one needs to know the effective Lagrangians (and the 
coupling constants appearing therein) at (a) the nucleon-nucleon-meson, (b) 
the resonance-nucleon-meson, and (c) the resonance- i^'''-hyperon vertices. 

K^ Y P 



N' 



JI,p,0,(X) 




Figure 1: Feynman diagrams 
for K^Y production in pp colli- 
sions. The diagram on the left 
shows the direct process while 
that on the right the exchange 
one. 



The parameters for A'^A'^ vertices are determined by fitting the A'^A^ elastic 
scattering T matrix with an effective A'^A^ interaction based on the vr, p, uj 
and a meson exchanges. The effective meson- A^A^ Lagrangians are 
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We have used the notations and conventions of Bjorken and Drell |32]. In 
Eq. (1) rn-Tv denotes the nucleon mass. Note that we have used a PV coupling 
for the NNtt vertex. Since we use these Lagrangians to directly model 
the T-matrix, we have also included a nucleon-nucleon-axial-vector-isovector 
vertex, with the effective Lagrangian given by 
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where A represents the axial-vector meson field. This term is introduced 
because in the limit of large axial meson mass (rriyi) it cures the unphysical 
behavior in the angular distribution of NN scattering caused by the contact 
term in the one-pion exchange amplitude [HH]) if gNN A is chosen to be 
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with very large (» m-Ar) rriA- fw appearing in Eq. (6) is related to qnNtt as 
We introduce, at each interaction vertex, the form factor 



Fr = {^^^^^],i = n,p,a,u;, (7) 



where qi and ?7i, are the four momentum and mass of the zth exchanged 
meson, respectively. The form factors suppress the contributions of high 
momenta and the parameter Aj, which governs the range of suppression, 
can be related to the hadron size. Since NN elastic scattering cross sections 
decrease gradually with the beam energy (beyond certain value), we take 
energy dependent meson-nucleon coupling constants of the following form 

g{y/s) = 9o(ixp{-iy/s), (8) 

in order to reproduce these data in the entire range of beam energies. The 
parameters, qq, A and H. were determined by fitting to the elastic proton- 
proton and proton-neutron scattering data at the beam energies in the range 
of 400 MeV to 4.0 GeV [SIIIHI- It may be noted that this procedure fixes 
also the signs of the effective Lagrangians [Eqs. (l)-(5)]. The values of 
various parameters are given in Table 1 of Ref . ^21, . The same parameters 
for these vertices were also used in calculations of the pion and the dilepton 
production in pp collisions. Thus we ensure that the NN elastic scattering 
channel remains the same in the description of various inelastic channels 
within this approach, as it should be. 

Below 2 GeV center of mass (cm.) energy, only three resonances, 
A^*(1650), iV*(1710), and iV*(1720), have significant decay branching ra- 
tios into KY channels. Therefore, we have considered only these three 
resonances in our calculations. The A^*(1700) resonance having very small 
(and uncertain) branching ratio for the decay to these channels, has been ex- 
cluded. Since all the three resonances can couple to the meson-nucleon chan- 
nel considered in the previous section, we require the effective Lagrangians 
for all the four resonance-nucleon-meson vertices corresponding to all the 
included resonances. Since the mass of the strange quark is much higher 
than that of the u— or d— quark, one does not expect the pion like strict 
chiral constraints for the case of other pseudoscalar mesons like r/ and K (to 
be called Q in the following). Thus, one has a choice of psuedoscalar (PS) or 
PV couplings for the NNC, and N'^.^NQ vertices (forms of the corresponding 
effective Lagrangians are given in Ref. !2lJ. The same holds also for the 
Ny^YK vertices. 



In principle, one can select a linear combination of both and fit the 
PS/PV ratio to the data. However, to minimize the number of parameters 
we choose either PS or PV coupling at a time. In the results shown below, 
we have used PS couplings for both N*Ntt and N*AK~^ vertices involving 
spin-1/2 resonances of even and odd parities. Calculations have also been 
performed with the corresponding PV couplings. The cross sections calcu- 
lated with this option for the N^i^YK vertex deviate very little from those 
obtained with the corresponding PS couplings. However, data shows a clear 
preference for the PS coupling at the N^.^Ntt vertices. 

The couplings constants for the vertices involving resonances are deter- 
mined from the experimentally observed quantities such as branching ratios 
for their decays to corresponding channels. It may however, be noted that 
such a procedure can not be used to determine the coupling constant for 
the A^*(1650)Si^ vertices, as the on-shell decays of this resonance to T,K 
channel are inhibited. Therefore, we have tried to determine this coupling 
constant by fitting the available data on the vr+p — > T,~^K~^, tt~p -^ T,^K^, 
and 7T~p -^ Ti~K~^ reactions in an effective Lagrangian coupled channels ap- 
proach inHEHl) where all the available data for the transitions from vrA^ to 
five meson-baryon final states, vrA^, ttttN, r]N, KA, and KY, are simultane- 
ously analyzed for center of mass energies ranging from threshold to 2 GeV. 
In this analysis all the baryonic resonances with spin < | up to excitation 
energies of 2 GeV are included as intermediate states. Since the resonances 
considered in this study have no known branching ratios for the decay into 
the A^a; channel, we determine the coupling constants for the N*Nuj vertices 
by the strict vector meson dominance (VMD) hypothesis ^^, which is based 
essentially on the assumption that the coupling of photons on hadrons takes 
place through a vector meson. 

It should be stressed that the branching ratios determine only the square 
of the corresponding coupling constants; thus their signs remain uncertain 
in this method. Predictions from independent calculations {e.g the quark 
model) can, however, be used to constrain these signs. The magnitude as 
well as signs of the coupling constants for the N*Ntt, N*AK, N*Np, and 
N* N (mr) s-wave vertices were determined by Feuster and Mosel jHl] and 
Manley and Saleski j^Z] in their analysis of the pion-nucleon data involving 
the final states vrA^, ttttN, r]N, and KA. Predictions for some of these quan- 
tities are also given in the constituent quark model calculations of Capstick 
and Roberts ^2Sj- Guided by the results of these studies, we have chosen the 
positive sign for the coupling constants for these vertices. Unfortunately, 
the quark model calculations for the N*Nlv vertices are still sparse and an 
unambiguous prediction for the signs of the corresponding coupling con- 



stants may not be possible at tliis stage [SHI- Nevertlieless, we liave cliosen 
a positive sign for the coupling constants for these vertices as well. 

The resonance properties used in the calculations of the decay widths are 
given in Table 1, where the resulting coupling constants and the adopted 
values of the cut-off parameters (A^ ) are also shown. It may be noted 
that we have fixed the latter to one value in order to minimize the number 
of free parameters. 

Table 1: Coupling constants and cut-off parameters for the A^*A^- meson and 
iV*-hyperon-meson vertices used in the calculations 



Resonance 


Decay channel 


Branching 
ratio (GeV) 


g-/A^ 


cut-off 

(GcV) 


N*{nio) 


A^vr 


0.150 


0.0863 


850.0 




Np 


0.150 


1.3653 


850.0 




Nuj 




0.1189 


850.0 




Na 


0.170 


0.0361 


850.0 




AK 


0.150 


2.9761 






SET 




4.4044 




A^*(1720) 


A^vr 


0.100 


0.0023 


850.0 




Np 


0.700 


90.637 


850.0 




Nuj 




22.810 


850.0 




Na 


0.120 


0.1926 


850.0 




AK 


0.080 


0.0817 






T.K 




0.2204 




Af*(1650) 


Ntt 


0.700 


0.0521 


850.0 




Np 


0.08 


0.5447 


850.0 




Nlo 




0.2582 


850.0 




Na 


0.025 


0.2882 


850.0 




AK 


0.070 


0.0485 






-EK 




0.0165 





After having established the effective Lagrangians, coupling constants 
and form of the propagators (which are given in Ref. j21|), it is straight 
forward to write down the amplitudes for various diagrams associated with 
the pp -^ pYK reactions which can be calculated numerically by following 
e.g. the techniques discussed in J2H]- The isospin part is treated separately. 
This gives rise to a constant factor for each graph, which is unity for the 



reaction under study. It should be noted that the signs of various ampUtudes 
are fixed by those of the effective Lagrangian densities, couphng constants 
and propagators as described above. These signs are not ahowed to change 
anywhere in the calculations. 

In the present form of our effective Lagrangian theory, the energy depen- 
dence of the cross section due to FSI is separated from that of the primary 
production amplitude and the total amplitude is written as, 

Afi = Mf,{pp^pYK+)-Tff, (9) 

where Mfi{pp — > pYK^) is the primary associated YK production ampli- 
tude, while Tff describes the re-scattering among the final particles which 
goes to unity in the limit of no FSI. The latter is taken to be the coherent 
sum of the two-body on-mass-shell elastic scattering amplitudes tj (with i 
going from 1 to 3), of the interacting particle pairs j — fe in the final channel. 
This type of approach has been used earlier to describe the pion [29 l lin | l^. 
r/- meson |421 1431 El] . h.K~^ j21j and (/>- meson j45j production in pp collisions. 
An assumption inherent in Eq. (9) is that the reaction takes place over 
a small region of space (which is fulfilled rather well in near threshold re- 
actions involving heavy mesons). Under this condition the amplitudes tj 
can be expressed in terms of the inverse of the Jost function |271 I29j which 
has been calculated by using a Coulomb modified effective range expansion 
of the phase-shift [IHI- The required effective range and scattering length 
parameters are given in Refs. |211l22j . 

3 Results and Discussions 

The total cross sections for the KK~^ and TPK^ reactions as a function of 
the excess energy are shown in Fig. 2. The calculations are the coherent 
sum of all resonance excitation and meson exchange processes as described 
earlier. In both cases a good agreement is obtained between theory and 
the data available from the COSY-11 collaboration. Keeping in mind the 
fact that all parameters of the model, except for those of N*Yp vertices 
and the FSI, were the same in the two calculations and that no parameter 
was freely varied, this agreement is quite satisfactory. It should be noted 
that we do not require to introduce arbitrary normalization constants to 
get the agreement between calculations and the data. We also show in this 
figure the results obtained without including the FSI effects (dashed line). 
It can be seen that the FSI effects are vital for a proper description of the 
experimental data in both the cases. 
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Figure 2: Comparison of the 
calculated and the experimental 
total cross sections for the pp — > 
ptsJi^ and pp -^ pTF'K^ reac- 
tions as a function of the excess 
energy. Results obtained with no 
FSI effects are shown by dashed 
lines. The experimental data are 
from Refs. [TT] . 
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In Fig. 3, we have investigated the role of various meson exchange processes 
in describing the total cross sections. The dashed, long-dashed, dashed- 
dotted, and solid with black square curves represent the contributions of vr, 
/9, Lo and a meson exchanges, respectively. The contribution of the heavy 
axial meson exchange is not shown in this figure as it is negligibly small. The 
coherent sum of all the meson-exchange processes is shown by the solid line. 
It is clear that the pion exchange graphs dominate the production process 
for both the reactions in the entire range of beam energies. Contributions of 
p and LO meson exchanges are almost insignificant. On the other hand, the 
a meson exchange, which models the correlated s— wave two-pion exchange 
process and provides about 2/3 of this exchange in the low energy NN 
interaction, plays a relatively more important role. This observation has 
also been made in case of the NN -^ NNtt reaction [I7 | UHl liOl [28| . 

The individual contributions of various nucleon resonances to the total 
cross sections of the two reactions are shown in Fig. 4. We note that in both 
the cases, the cross section is dominated by the A^*(1650) resonance excita- 
tion for e < 30 MeV. Since A'^*(1650) is the lowest energy baryonic resonance 
having branching ratios for the decay to YK~^ channels, its dominance in 
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Figure 3: Contributions of vari- 
ous meson exchange processes to 
the total cross section for the 
PP ^ pK^A and pp -^ pK^T,'^ 
reactions. The dashed, long- 
dashed, dashcd-dotted and solid 
with black squares curves repre- 
sent the contributions of tt, p, to 
and (7 meson exchanges, respec- 
tively. Their coherent sums are 
shown by the solid lines . 
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these reactions at beam energies very close the kaon production threshold 
is to be expected. In the near threshold region the relative dominance of 
various resonances is determined by the dynamics of the reaction where the 
difference of about 60 MeV in excitation energies of A^*(1650) and A^*(1710) 
resonances plays a crucial role. 

However, for e values between 30 - 60 MeV, while the pp -^ pK~^A reac- 
tion continues to be dominated by A^*(1650) excitation, the pp -^ pK'^YP 
reaction gets significant contributions also from A^*(1710) and A^*(1720) 
resonances. This difference in the role of the three resonances in the two 
cases can be understood in the following way. For a resonances to con- 
tribute significantly, we should have my -|- ttik + e > mn + Tr/2, where 
m/j and Tn are the mass and width of the resonance, respectively. There- 
fore, in the region of excess energies > Q[= (m^ -|- Tii/2) — {my + m,/^)], 
the particular resonance R should contribute significantly. The values of Q 
for the pp -^ pK~^A reaction, are 115 MeV, 150 MeV, and 185 MeV, for 
the iV*(1650), iV*(1710), and iV*(1720) resonances, respectively. On the 
other hand, for the pp -^ pK^YP case, they are 36 MeV, 72 MeV and 105 
MeV, respectively for these three resonances. Therefore, contributions of 
the A^*(1710) and A^*(1720) resonance excitations are negligibly small for 
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the A'^A production in the entire energy range of the COSY- 11 data (i.e., 
for e < 60. MeV) while they are significant for the K^TP case for e > 30 
MeV. It would be helpful to have data on the invariant mass spectrum at 
these excess energies in order to confirm these theoretical observations. 
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Figure 4: Contributions of N*(1650) (dashed line), N*(1710) (full line with black 
squares) and N*(1720) (dashed-dotted line) baryonic resonances to the total cross 
section for the two reactions studied in Fig. 2. Their coherent sum is shown by the 
solid line. 

In Fig. 5, we compare our calculations with the data for the ratio i? as a 
function of e. We have shown here the results for excess energies up to 60 
MeV, where the COSY-11 data are available. It is clear that our calcula- 
tions are able to describe well the trend of the fall-off of i? as a function 
of the excess energy. It should be noted that FSI effects account for about 
60% — 80% of the observed ratio for e < 30 MeV and about 50% of it beyond 
this energy. Therefore, not all of the observed value of R at these beam en- 
ergies can be accounted for by the FSI effects, which is in agreement with 
the observation made in ^Hj. It should again be emphasized that without 
considering the contributions of the A^*(1650) resonance for the Y?K^ re- 
actions, the calculated ratio would be at least an order of magnitude larger. 
Therefore, these data are indeed sensitive to the details of the reaction mech- 
anism. At higher beam energies (e > 300 MeV), values of R obtained with 
and without FSI effects are almost identical. In this region the reaction 
mechanism is different; here all the three resonances contribute in one way 



11 



or the other, their interference terms are significant |2J, and FSI related 
effects are unimportant. 
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Figure 5: Ratio of the total cross sections for pp — > pAK'^ and pp -^ pTPK^ 
reaction as a function of the excess energy. The solid and dashed lines show the 
results of our calculations with and without FSI effects, respectively. The data are 
from 17 . 



4 Summary and conclusions 

In summary, we have studied the pp -^ pAK^ and pp — > pT,^K~^ reactions 
within an effective Lagrangian model. Most of the parameters of the model 
are fixed by fitting to the elastic A^A^ T-matrix which restricts the freedom 
of varying them freely in order to fit the data. The reactions proceed via 
excitation of the N*{1650), A^*(1710), and iV*(1720) intermediate baryonic 
resonant states. An important result of our study is that the A^*(1650) 
resonant state contributes predominantly to both these reactions at near 
threshold beam energies. Therefore, these reactions in this energy regime, 
provide an ideal means of investigating the properties of this 5ii baryonic 
resonance. To the extent that the final state interaction effects in the exit 
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channel can be accounted for by the Watson-Migdal theory, our model is 
able to explain the experimentally observed large ratio of the total cross 
sections of the two reactions in the near threshold region. 
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